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We perform an exact enumeration study of polymers formed from a (quenched) random sequence 
of charged monomers ±go, restricted to a 2-dimensionaI square lattice. Monomers interact via a 
logarithmic (Coulomb) interaction. We study the ground state properties of the polymers as a 
function of their excess charge Q for all possible charge sequences up to a polymer length A'^ = 18. 
We find that the ground state of the neutral ensemble is compact and its energy extensive and 
self- averaging. The addition of small excess charge causes an expansion of the ground state with the 
monomer density depending only on Q. In an annealed ensemble the ground state is fully stretched 
for any excess charge Q > 0. 
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I. INTRODUCTION 



Polymer systems have been subject to extensive study 
for several decades The interest in such macro- 

molecules arises from many sources, where one impor- 
tant reason is the natural occurrence of polymers in bio- 
logical systems, for example proteins. Random polymers 
are also good examples of disordered systems, sharing 
common features with other physical systems exhibit- 
ing randomness or frustration An important class 
of polymers are polyampholytes (PA), where the chain 
is composed of positive and negative charged monomers 
1^. Certain properties of proteins, some of whose con- 
stituents, the amino acids, may be charged may be 
understood by studying PAs. 

Much attention has been devoted to studying the 
ground state properties of polymers (or the native state 
of proteins) and understanding which features determine 
their shape ^,||. In PAs, the long-range Coulomb in- 
teraction between the monomers may, in the absence 
of screening, control the size and shape of the poly- 
mer, dominating over other short-range interactions. For 
d > 4, the behavior of the PA is similar to that of ran- 
dom chains with short-range interactions, the Coulomb 
interaction becoming an irrelevant perturbation For 
d < 4, the effect of the Coulomb interaction is compli- 
cated. Higgs and Joanny and Wittemer et al. [ p^ , 
expanding on arguments previously made by Edwards et 
al. [|ri|, pointed out, on the basis of the Debye-Hiickel 
theory |l2|, that charge fluctuations along the chain in- 
duce a net attractive interaction between segments of the 
chain, leading to a collapse of the PA to a compact shape. 
The resulting picture is of a globule consisting of densely 
packed blobs. Each blob is a dilute region of a weakly 
perturbed self-avoiding walk (SAW) of the size of the 
Debye-Hiickel screening length. Lowering the tempera- 
ture, decreases the size of the blobs. This behavior is not 
a phase transition but rather a gradual increase of den- 
sity of monomers with decreasing temperature. Charge 
correlations along the polymer may change the nature of 
the collapse llQ] : a polymer with a quenched sequence of 



alternating charges will undergo a ^-transition |13[, like 
polymers with short-range attraction interaction |14[. 

The application of the Debye-Hiickel theory requires 
the strict neutrality of the system Q. If we assign to 
each monomer along the chain a charge qi = iqo, then 
an excess charge Q = 9* be defined. A ran- 

dom PA will not necessarily be neutral but may carry an 
excess charge Q, typically on the order of qoVN. Gutin 
and Shakhnovich jl^ , and Dobrynin and Rubinstein |l^] 
treated the problem of a charged PA. They found an ex- 
pansion of the ground state due to excess charge and 
describe it in the terms of an elongated globule with a 
Q-dependent aspect ratio. Observing that for neutral 
quenches {Q = 0) the ground state is compact, Kantor 
and Kardar ||l^,|l^ described the ground state proper- 
ties of a charged PA in an analogy to a drop of a fluid 
with a charge distributed in it: The energy of a neutral 
drop may be described by an extensive form with a sur- 
face tension accounting for the compact spherical shape. 
The excess charge Q is the important parameter deter- 
mining the properties of the drop (or compact PA). A 
phenomenological description of the energy is given by 



E ^ ^ecV + jS + Q^F{R) 



(1) 



where the first term is the extensive contribution to the 
energy (proportional to the volume V), the second term 
is a correction proportional to the surface area S, while 
the third term is the Coulomb interaction, where F{R) 
represents the dependence of the Coulomb potential on 
the linear size of the drop R and is given (up to dimen- 
sionless prefactors of order unity) by: 



^ ' 1 -Ini? ford = 2. 



(2) 



Increasing the charge Q creates an outward pressure that 
seeks to expand the drop, competing with the surface 
tension. At a certain excess charge, called the Rayleigh 
charge Qn, the pressure difference between the inside and 
outside of the drop vanishes and the drop becomes locally 
unstable to elongation. In a d-dimensional PA this hap- 



pens at Qjj 



18(1. However, even at a lower 
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excess charge it might be energeticahy favorable to split 
a drop into two distant droplets each carrying half the 
excess charge leading to a global instability. The critical 
charge Qc when this occurs is found to be — aQj^, 
where the proportionality constant depends on d. This is 
the excess charge above which the PA should expand. In 
3D, both the typical excess charge and the critical charge 
scale as V^, hence, a typical PA will be stretched. Nu- 
merical studies using Monte-Carlo methods and ex- 
act enumeration techniques support these qualitative 
predictions in 3D. Once the critical charge is exceeded, 
the PA expands, but it does not have the freedom of the 
drop to disintegrate, since the polymer is a connected 
object. It can, however, be approximately described as 
a necklace jl^ of weakly charged globules connected by 
highly charged strands of a PA. 

While the 3D case is important from the practical point 
of view, certain aspects cause conceptual difficulties (and 
consequently difficulties in the interpretation of Monte- 
Carlo results): In 3D, both Q^j and the critical charge 
which determines the high-temperature behavior of a PA 
and the typical fluctuation of the excess charge Q be- 
tween the quenches, have the same scaling - all three 
quantities are proportional to ^/N. It is, therefore, bene- 
ficial to study PAs in other space dimensions d. We study 
a model of a PA, a 2-dimensional (2D) randomly charged 
SAW on a 2D square lattice. Although this model signifi- 
cantly simplifies the complexity of real physical systems, 
we hope to gain some physical insight and understand 
the relevant mechanisms governing the behavior of the 
system. In Sec. || we elaborate on the specific model we 
study. 

The logarithmic behavior of a 2D Coulomb potential 
creates difficulties which are absent in d > 2: It is even 
not self-evident that the low temperature configurations 
have an extensive energy and surface tension, i.e. that 
the description of Eq. (|^) is applicable in 2D. In Sec. |l| 
we will show that the extensive description of the energy 
applies to our 2D model for Q ~ 0. In Sections ^,[v| we 
study the effect of excess charge on 2D PAs. Although 
the short chains studied prevent the verification of the 
necklace description, we see that any deviation from neu- 
trality leads to an expansion of the ground state. We find 
that the system is much more susceptible to excess charge 
in 2D than in 3D. 



II. THE MODEL AND METHOD 

We model the polymer as a SAW on a 2D square lat- 
tice of spacing a. Monomer i along the chain has a charge 
Qi which is randomly assigned a value ±go- The charges 
interact via a 2D Coulomb potential 

E^-y2q,q,\n^ , (3) 



where is the distance between the zth and jth 
monomers and Tq is an arbitrary constant setting the 
reference point of the energy. The study of a 2D system 
allows us to extend our study to longer polymers than 
was possible in 3D |]l8|. In addition, for short 2D chains, 
the surface to volume ratio, responsible for many of the 
properties of the poly mer, is similar to that of l ong , re- 
alistic 3D chains |2^. Everywhere (except Sec. we 
consider a fixed charge sequence, i.e. quenched disorder. 
Note that the energy (||) is symmetric under charge con- 
jugation, so we study only overall neutral and positively 
charged quenches. 

In this study we apply the complete enumeration 
method: We create all possible spatial configurations of 
the chain and calculate exact quantities for each config- 
uration. This allows us to find the ground state of the 
system and, when studying higher temperatures (results 
we do not report here) , to compute exact thermodynamic 
averages. While the number of SAWs grows exponen- 
tially with N Q limiting us to relatively short chains, 
this is the only method for studying exact ground state 
properties, since Monte Carlo methods fail to equilibrate 
at low temperatures. Up to iV = 18 we enumerated 
all possible charge sequences. For a PA 18 monomers 
long, after taking into account basic symmetries, there 
are 5,808,335 spatial configurations for each of the 77,819 
possible quenches (unrelated by symmetry). Enumerat- 
ing all configurations for all the quenches of = 18 
required approximately 20 days of CPU time on a Sil- 
icon Graphics RIOOOO workstation. Up to a length of 
26 monomers we enumerated only partial samples of all 
possible quenches (1000, 500, 100 and 50 random neutral 
quenches for lengths N — 20, 22, 24 and 26 respectively). 
A 26 monomer long chain has 15,435,169,364 possible 
spatial configurations, unrelated by symmetry. Similar 
model and method were employed by Kantor and Kar- 
dar to study 3D systems up to A^ = 13. 

There is an important point to note regarding the en- 
ergy of a 2D system with a fixed number of charges. We 
may express the distances as = ar[j , where r[j is a di- 
mensionless distance (in lattice constants), and re- write 
Eq. (^) in the form: 

E^-Y.^^q, In -\{Q^- qlN) In - . (4) 

The first term depends on the specific spatial conforma- 
tion. The second term is a constant reference point, in- 
dependent of the spatial configuration or charge sequence 
along the polymer, and is a function of the excess charge 
only. Here on, we take tq — a, so that the second term 
vanishes and calculate only the first term as the "en- 
ergy." Within an ensemble of the same excess charge, 
this has no effect since we are interested only in energy 
differences. However, it must be kept in mind that for 
different excess charges and lengths, the choice of rg may 
set different reference points, so energies cannot be com- 
pared. 
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III. GROUND STATE SPECTRUM OF NEUTRAL 
POLYAMPHOLYTES 



Following the reasoning given in Sec. | and empirical 
observations we expect the ground state of neutral PAs 
to be compact. It would be natural to consider an exten- 
sive form of the energy of the ground state with a surface 
correction. Phcnomenologically, this can described as: 



(5) 



where represents an average over quenches. The sec- 
ond term in represents the 2D surface correction. Un- 
like short-range interacting systems, it is not obvious that 
the 2D Coulomb interaction produces such a form of E. 
For instance, the energy of a randomly selected spatial 
configuration of a neutral sequence is expected to be of 
order Nh\N, i.e. it grows faster than N. 

A convenient comparison point for the ground state en- 
ergies is provided by a checkerboard configuration ( "salt 
crystal") of charges, which is the ground state of a 2D 
Coulomb gas on a lattice |2^]. It can be shown that for 
such an arrangement of charges the energy is extensive. 
Dots in Fig. ^ depict the energy per monomer for finite 
"salt crystal" configurations forming a square consisting 
oi N = An? {n =2, 3, 4, ...) charges. The energy of such 
a crystal is given by 



E= -OMq^N 



0.26ql 



N 



(6) 



Not surprisingly, in PAs the ground state of quenches 
which are an alternating sequence of positive and nega- 
tive charges (open circles in Fig. |l]), follow the same law. 
However a typical quench cannot assume a salt-crystal- 
like ground state. Results of evaluating all ground state 
energies for the different quenches are depicted in Fig. 
|l|. For each length, we display the mean and standard 
deviation of the ground state energy taken over the dif- 
ferent quenches. The mean energies may be fitted to an 
extensive form: 



E = (-0.133 ± 0.005)q^N + (0.30 ± 0.02)q^VN . (7) 



Results must be considered with caution due to the 
rather short chains studied. A definite exclusion of possi- 
ble logarithmic corrections would require the considera- 
tion of chain lengths larger by at least one order of mag- 
nitude. Nevertheless, since the energy of a salt crystal 
bounds from below the ground state energies, we believe 
that the extensive form is valid. 
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FIG. 1. Ground state energies per monomer ||28|| for neutral 
sequences of lengths 6-26 monomers. For each length, boxes rep- 
resent the mean value over all quenches studied, while error bars 
represent the standard deviation. The results for the alternating 
sign quench are represented by open circles. The energy of a "salt 
crystal" plane is represented by dots. 
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Ground State Energy Per Monomer 
FIG. 2. A histogram of the ground state energy per monomer 
psi for all quenches of lengths 16 and 18 monomers and a partial 
sample of quenches of length 20 monomers. Crosses indicate mean 
and standard deviation limits. Bin size is O.Olgp. 



The distribution of the ground state energies per 
monomer, for the ensemble of quenches, is relatively 
broad even for the largest Ns considered: Fig. |^ depicts 
the distribution of the ground state energies of =16 
and 18 (all possible quenches), and 20 monomers (a par- 
tial sample only). Nevertheless, the standard deviation a 
of the energy per monomer gradually decreases with in- 
creasing A^. In systems with short-range interactions the 
energy of a large random system may deviate from the 
ensemble-averaged mean by an amount \/N. This is, 
usually, a consequence of the fact that the total energy of 
the system is (approximately) a sum of the energies of the 
subsystems: E.g., the total energy is a sum of energies 
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obtained by dividing the system into two halves. This 
property is referred to as self-averaging. Fig. ^ demon- 
strates the linear dependence of the standard deviation of 
the energy per monomer cr on 1/ V^, which means that 
the fluctuations of the total energy are ~ \/N . Thus, we 
may treat our system as self-averaging. However, if we 
split a large 2D neutral PA into two equal parts, usually 
none of the halves will be neutral. On the contrary, typ- 
ically each half will be charged by an equal and opposite 
charge Q' ~ Vn , and, thus, have electrostatic energy of 
order — Q'^lni?' ^ —N\nN. Moreover, as will be shown 
in the next section, such charged sub-chains will be ex- 
panded, i.e. the spatial conformations of the sub-chains 
will not resemble the configuration of the entire PA. Nev- 
ertheless, we may use a slightly modified concept: In a 
very long neutral PA, one can, usually, find several special 
points, such that splitting the the PA at those points will 
divide it into several large newira/ sub-systems (see, e.g., 
Rcf . ||2^ ) . The total energy of a PA can now be regarded 
as the sum of the energies of these separate segments. 
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FIG. 3. Standard deviation of the ground state energy per 
monomer a , as a function of 1/y/N, for neutral quenches of 
lengths 6, 8,. ..,26 monomers. Error bars denote 98% confidence 
levels when the standard deviation was estimated for partial sam- 
ples, in the four leftmost cases. 



IV. EXCESS CHARGE EFFECTS 

Non-vanishing excess charge Q has a profound effect 
on the energies and spatial conformations of PAs. Fig. ^ 
depicts the ground state energies of all quenches of PAs of 

= 18. Every point in this figure corresponds to a differ- 
ent quench. The ground state energies split into bands, 
each of which corresponds to a different Q. The dis- 
tance between the band of neutral quenches and the rest 
of the bands increases approximately as due to the 
appearance of a long-range electrostatic term — Q^lni?, 
where the length R characterizes the spatial extent of 
the PA in the ground state. (The length R depends on 
Q, however, Ini? has a minor effect on the leading 



dependence.) One should keep in mind that, in general, 
the absolute values of the energy difference between the 
bands depends on the choice of the reference point rp in 
Eq. (^. We note that bands remain relatively narrow, 
i.e. the ground state energy weakly depends on the de- 
tails of the specific quench. Only when Q approaches 
q^N , the bands broaden, i.e. the details of the quench 
have a more significant inffuence on the overall energy. 
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Arbitrary Order of Quenches (xlO*) 

Ground state energies for all quenches of length 18 
monomers [ p^ . Abscissa depicts an arbitrary order of quenches. 
Annotations denote the excess charge of the bands the energies 
group up into. 
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Excess Charge - Q 
FIG. 5. R^/N vs. Q [ p8| . Curves represent different lengths 
(A'^ =6, 7,. ..,17,18 from bottom to top). Note the collapse of the 
curves at low Q and gradual deviation as Q is increased. 

Excess charge is an important parameter determining 
the shape of the ground state as well. Specifically, it 
determines whether the ground state is compact or ex- 
panded. We measure the typical size of a configuration 
by its root-mean-square size called the radius of gyration 
Rg. A compact shape is defined by Rg ^ N^/'''. In 2D 
this means that Rg/N is independent of N. In Fig. || we 
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plot Rg/N as a function of Q for the different values of 
N. Each point is an average over all possible quenches of 
the specific N and Q. The most striking feature of Fig. 
H is that for small Q all the curves, for different values 
of N, collapse onto the same curve. Thus, even for very 
large N a minute charge Q can significantly increase Rg. 
This means that among the ensemble of random quenches 
in the N ^ oo limit a vanishing portion remains com- 
pact. Even if the ensemble is biased towards neutrality, 
i-e. Qtypicai ~ N^, where 1 a; > 0, the resulting con- 
figurations will be expanded. Only exceptional quenches, 
where the excess charge is very small even for large N, 
will remain compact. This result may also be incorpo- 
rated within the charged drop model for a PA. In 2D 
Qr ^ iV^/'*, but the logarithmic potential creates a cer- 
tain difference with respect to higher dimensions: Be- 
cause the potential has no finite asymptotic value, the 
charged drop will always find it favorable to disintegrate 
into distant droplets with an infinite energy gain. Hence, 
in a continuum model Qc — 0, and we except an expan- 
sion in the size of the ground state of the PA for any 
excess charge. This is very different from 3D systems 
where Qc (x Qr \/N so a finite portion of 
all random quenches remains compact. Once the critical 
charge is exceeded, the ground state will begin expanding 
but does not immediately become fully stretched. Nu- 
merical results show that in 2D for Q ^ A^"-^ the result- 
ing configuration is completely stretched, i.e. Rg ^ N . 
This implies that if the charge is a finite fraction of the 
size, i.e. Q ^ N, the ground state is fully stretched. 



V. THE ANNEALED ENSEMBLE 



Thus far we have considered quenched sequences in 
which charges are fixed in position along the chain. We 
now remove this restriction, allowing charges to exchange 
positions along the chain, maintaining the overall excess 
charge and number. (Two charges cannot be located 
at the same site, thus creating a double or vanishing 
charge.) We shall refer to this as an annealed ensemble of 
PAs. The ground state configurations for the ensembles 
of A^ = 16 and different values of Q are shown in Fig. ||. 

For the neutral ensemble {Q = 0) the ground state is 
naturally a sequence of charges alternating in sign that 
form a "salt crystal" configuration, in the same way that 
free charges form such a crystal This spatial ar- 

rangement of the charges can be obtained by many differ- 
ent spatial configurations of the chain, making the ground 
state highly degenerate. For any non-zero excess charge 
we observe a sharp transition and the ground state has 
a fully stretched rod-like shape. For the sizes consid- 
ered here, the ground state configuration and charge se- 
quences were unique. 
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FIG. 6. Ground state configuration for annealed sequences of 
length 16, with an overall excess charge restriction Q. Opposite 
charge types are denoted by open and full circles. 

Previously, we have seen that for a finite Q, the ground 
state configuration of the quenched ensemble expands. 
In the annealed case the expansion is much more abrupt: 
The PA is completely stretched for any minute Q. It was 
suggested that an annealed PA in 3D, with a large 
enough excess charge, might find it favorable to expel 
the excess charge in the form of highly charged fingers, 
while most of the monomers remain in a globule. In 2D 
the behavior seems to be quite different: Small excess 
charge suffices to stretch the PA into a straight line. The 
positions of the charges on that line are determined as 
follows: in general the charges attempt to maintain an 
alternating sequence; this not strictly possible because of 
nonvanishing Q. The excess charge is spread out almost 
uniformly (up to logarithmic corrections) along the PA. 

These results should be regarded with caution. Con- 
trary to the quenched ensemble, where the ground state is 
averaged over many quenches, thus smoothing out lattice 
effects, the ground state of the annealed ensemble corre- 
sponds to a single configuration of a single sequence. The 
chains studied might not be long enough to overcome this 
problem. 



VI. DISCUSSION 

We studied the ground state properties of a model 
2D PA. We found that neutral quenches have a com- 
pact, dense ground state, whose energy can be described 
as an extensive energy term with a surface correction. 
The ground state energies of the different quenches may 
be described by energy bands, each band correspond- 
ing to a certain excess charge. The specific sequence 
of charges is of a small influence, creating small fluc- 
tuations about the mean within the bands (for neutral 
quenches the ground state energy is even self-averaging 
over all possible quenches). When the charge sequence is 
annealed (not quenched) the shape of the ground state 
is much more sensitive to excess charge, i.e. the PA is 
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fully stretched for any Q ^ Q. The results are differ- 
ent than those obtained for a similar 3D model [^ . 
The nature of the 2D Coulomb interaction accounts for 
this difference. Both cases, however, match the expecta- 
tions derived from an analogy with a charged drop. The 
main drawback of the method we applied throughout this 
study is the limited size of the systems we can fully in- 
vestigate. Although we surpassed the maximal length 
studied in 3D, some results may be suspected as arising 
from finite size and lattice effects. A major increase in 
computer capabilities is required to carry this method to 
larger systems. 

Several authors jl^,^ suggested that, at low temper- 
atures, a PA may undergo a freezing transition as do 
some models of polymers with a short-range interaction 
|p|, p4|j2^ and other disordered systems. Should this be 
the case, then the ground state properties do not provide 
the full description of the zero-temperature behavior of 
the system. Our investigation of such a transition was 
not conclusive, but indicates its absence, at least in the 
sense that the requirements of the random energy model 
p6t are violated. This is in agreement with predictions 
of Pande et al. that long-range interacting systems 
would not undergo a glass transition. 
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